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ABSTRACT 


The  three  elastic  constants  of  tantalum,  tungsten  and  molybdenum 
have  been  measured  in  the  temperature  range  between  4.2  and  300°K.  The 
directly  measured  quantities  and  ^(C-q-C^)  f°r  molybdenum  and  the 
computed  bulk  moduli  of  all  three  "elements  show  an  anomalous  temperature 
dependence.  These  data  along  with  published  results  for  other  metals 
generally  show  a  difference  in  temperature  dependence  of  the  bulk  modulus 
between  3d  and  4d-5d  metals.  Overlap  contributions  to  the  elastic  con¬ 
stants  are  found  to  favor  an  antiferromagnetic  ordering.  The  Debye 
characteristic  temperatures  calculated  from  elastic  constants  are  found 
to  be  in  substantial  agreement  with  results  obtained  from  calorimetry. 
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INTRODUCTION 


The  elastic  constants  of  a  material  determine  the  linear  coefficients 
relating  stress  and  strain  in  that  material.  In  themselves  they  constitute 
an  important  equilibrium  property.  Because  the  equilibrium  properties  of 
a  system  are  related  through  considerations  of  the  system's  total  internal 
energy,  knowledge  of  the  elastic  constants  is  important  to  a  complete  un¬ 
derstanding  of  thermal,  electrical  and  magnetic  behavior. 

The  elastic  constants  of  a  solid  are  a  measure  of  the  curvature  of 
the  cohesive  energy  and  their  measurement  yields  information  about  the 
interatomic  forces.  In  particular,  low  temperature  observations  are  of 
importance  for  comparison  with  the  results  of  theoretical  calculations. 

Also,  the  low  temperature  elastic  constants  may  be  used  to  predict  the 
lattice  specific  heat,  which  may  then  be  compared  with  calorimetric  results. 

The  transition  metals  have  always  been  of  great  interest  because  of 
their  incomplete  d  shells.  The  Group  VB  and  Group  VIB  metals,  vanadium, 
chromium,  niobium,  molybdenum,  tantalum,  and  tungsten,  have  the  particular¬ 
ly  interesting  features  of  exhibiting  a  body-centered  cubic  structure  and 
very  high  cohesive  energies.  Until  recently  very  few  measurements  have 
been  made  on  these  metals,  principally  because  of  metallurgical  difficulties 
in  obtaining  pure  single  crystals. 
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EXPERIMENTAL  PROCEDURE 


The  adiabatic  elastic  constants  of  tantalum,  tungsten,  and  molybdenum 
have  been  obtained  from  room  temperature  down  to  liquid  helium  temperatures 
by  measuring  the  transit  times  of  elastic  waves  in  single  crystals.  Un- 

rectified  waveforms  were  observed  using  an  ultrasonic  apparatus  similar 

1  2 

to  that  described  previously.*  A  preliminary  report  on  the  elastic  con¬ 
stants  of  tantalum  pointing  out  anomalies  in  the  observed  bulk  modulus 
3 

has  already  been  given. 

The  tantalum  specimen  was  purchased  from  Materials  Research  Corpo¬ 
ration  as  a  quarter  inch  diameter  rod  from  which  a  2.0643  cm  (at  20°C) 
specimen  was  cut  with  accurately  parallel  faces.  Because  of  ifs  small 
diameter,  no  attempt  was  made  to  orient  this  crystal.  Instead  the  orien¬ 
tation  was  taken  "as  is".  This  orientation  was  determined  by  Laue  back- 
reflection  x-ray  methods  to  deviate  from  £lloJ  by  about  8°,  the  asso¬ 
ciated  polar  angles  being  9  =  84.9°  and  0  =  39.7°»  Transit  time  measure¬ 
ments  were  made  on  the  tantalum  specimen  both  before  and  after  annealing 
at  1200°C  for  five  hours  in  a  vacuum  of  at  least  10  ^  mm  Hg.  Annealing 
did  not  affect  the  observed  transit  times.  No  change  in  velocity  or  at¬ 
tenuation  of  the  elastic  waves  was  observed  in  passing  through  the  super¬ 
conducting  transition  temperature  for  tantalum.  The  transducer  binder 
used  with  the  tantalum  sample  was  glycerin  below  0°C  and  salol  at  and 
above  0°C. 

Temperatures  on  all  samples  were  controlled  to  within  0.3°C  during 
transit  time  measurements  using  a  Brown  Electronik  recorder.  A  copper- 
constantan  thermocouple  was  used  for  all  temperature  measurements  above 
liquid  air  temperatures.  For  temperatures  from  liquid  air  to  liquid 
helium  temperatures  a  copper  resistance  thermometer  was  used  when  making 
the  measurements  on  tantalum;  a  copper-gold  2.1  a/o  cobalt  alloy  ther¬ 
mocouple  was  used  for  the  tungsten  and  molybdenum  measurements  in  the 
same  low  temperature  range. 

The  tungsten  single  crystal  was  purchased  from  the  Linde  Company 
as  a  1/2"  diameter  rod  that  had  been  grown  with  its  longitudinal  axis 
2°  off  [j.loJ  •  From  this  a  specimen  was  prepared  that  had  two  parallel 
faces  accurately  perpendicular  to  [lio]  an(*  2.4866  cm  apart  (at  0°C)  . 
Nonaq  stopcock  grease  was  used  as  the  binder  for  tungsten  at  all  tem¬ 
peratures  . 
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The  molybdenum  single  crystal  was  obtained  from  the  Wah  Chang  Corpo¬ 
ration  as  a  rough  surfaced  bar  averaging  about  one  inch  in  diameter.  Using 
standard  Laue  back-reflection  techniques  a  QllcT]  direction  was  located 
in  this  specimen.  Two  parallel  faces  were  then  cut  perpendicular  to  this 
direction,  the  mean  distance  between  faces  being  1.8322  cm  at  0°C.  The 
maximum  deviation  from  parallelism  of  these  surfaces  in  the  molybdenum 
sample,  as  was  also  true  for  the  tantalum  and  tungsten  specimens,  was 
0.0002  cm.  The  transducer  binder  used  for  molybdenum  was  Nonaq  stopcock 
grease  for  the  readings  at  and  just  below  room  temperature  and  Dow  Corning 
High  Vacuum  Grease  from  below  these  temperatures  down  to  helium  temperatures. 

In  calculating  the  elastic  constants  of  these  metals  the  thermal  ex- 
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pansion  data  of  Nix  and  MacNair  were  used  down  to  the  limits  of  their  va¬ 
lidity  (the  liquid  air  temperature  region) .  Extrapolations  from  these  data 
following  Grueneisen  plots  were  used  down  to  4.2°K.  The  lattice  constant 
used  for  tantalum  was  3 „ 3058$  at  25°C,  as  given  in  a  survey  by  Swanson  and 

Tatge.  From  this  a  density  of  16.6324  grams  per  cubic  centimeter  at  0°C 

£ 

was  calculated.  For  tungsten  the  lattice  constant  given  in  Pearson  of 
3.1650&  at  25°C  was  used.  The  0°C  density  of  tungsten,  using  this  value 
of  lattice  parameter,  is  19.2710  grams  per  cubic  centimeter.  Pearson  was 
also  the  source  for  the  lattice  constant  of  molybdenum,  3.1468^  at  20 °C. 

The  0°K  density  of  molybdenum  calculated  from  this  value  is  10.2549  grams 
per  cubic  centimeter. 

Results  of  impurity  analyses  conducted  on  these  three  metals  are  con¬ 
tained  in  Table  I.  The  analyses  of  tantalum  and  molybdenum  were  conducted 
by  the  firm  of  Kennard  and  Drake  of  Los  Angeles,  using  test  material  that 
had  been  taken  from  the  original  specimens.  The  tungsten  analysis,  pro¬ 
vided  by  the  Linde  Company,  is  an  average  analysis  for  material  produced 
at  the  same  time  as  the  sample  used.  Table  I  shows  a  low  metal  impurity 
content  for  all  three  samples.  The  interstitial  gas  content  for  tungsten 
is  low,  larger  for  tantalum  and  quite  substantial  for  molybdenum. 
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TABLE  1=  Chemical  and  vacuum  fusion  analyses  of  crystal  specimens 
(ppm  by  weight) 


TANTALUM 

TUNGSTEN 

MOLYBDENUM 

Fe 

20 

C 

10 

Fe 

50 

Cu 

3 

S 

10 

Si 

77 

°2 

222 

A1 

<100 

°2 

930 

N2 

124 

Be 

<100 

n2 

140 

H2 

9 

Cd 

<100 

H2 

10 

Mg 

<1 

Si 

<100 

°2 

25 

N2 

<20 

H2 

1 

All  other  impurities  nil 
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EXPERIMENTAL  RESULTS 


For  the  tantalum  specimen  the  measured  transit  times  were  used  to 
2 

calculate  the  quantity  ^>V  for  each  of  the  three  modes  being  propagated 


in  a  direction  some  8°  off 
^  „2 


Smooth  lines  were  drawn  through 


these  raw  data  of  pV  versus  temperature,  yielding  the  curves  of  Fig.  1. 
Values  were  taken  from  these  curves  at  regular  increments  of  temperature 
and  a  successive  approximation  method  due  to  Neighbours^  was  applied  to 
derive  the  values  of  the  elastic  constants  listed  in  Table  II. 

Since  the  direction  of  elastic  wave  propagation  in  the  tungsten  and 
molybdenum  samples  was  directly  along  [lioj  >  the  values  of  calcu¬ 

lated  from  the  transit  time  measurements  are  equal  to  certain  linear  com¬ 
binations  of  the  three  adiabatic  elastic  constants.  Accordingly,  Fig.  2 

displays  the  measured  variation  with  temperature  for  the  longitudinal  mode 
2  1 

^>V  =  +  ^44  tungsten.  Fig.  3  displays  the  measured 

variation  with  temperature  of  the  ^>V^'s  associated  with  the  two  trans¬ 
verse  modes  for  tungsten.  These  equal  Y^ll"^12^  ^°r  t*ie  mo<*e  polarized 
perpendicular  to  the  z-axis  and  for  the  mode  polarized  parallel  to 

the  z-axis.  Smooth  curves  have  been  drawn  through  the  observed  data 
points  of  FigSo  2  and  3  and  a  compilation  is  presented  in  Table  III  of 
the  adiabatic  elastic  constants  of  tungsten  for  values  taken  at  regular 
intervals  of  temperature  from  these  smooth  curves « 

Correspondingly  Figs.  4  and  5  display  the  measured  longitudinal  and 
2 

transverse  pV  for  molybdenum.  Table  IV  shows  the  adiabatic  elastic 
constants  of  molybdenum  derived  from  the  curves  depicted  in  Figs.  4  and  5. 

Taking  into  account  the  various  sources  of  error,  we  estimate  the 

2 

absolute  accuracy  of  a  single  observation  of  ^>V  to  be  +  0.5  percent, 
although  the  internal  consistency  of  the  data  is  much  better  than  this. 
Since  the  elastic  wave  propagation  in  the  tungsten  and  molybdenum  samples 
was  directly  along  a  principal  direction,  the  above  estimate  of  +  0.5 
percent  should  be  valid  as  an  estimate  of  error  for  the  final  reported 
values  of  the  elastic  constants.,  In  the  case  of  tantalum,  however,  the 
direction  of  propagation  was  some  8°  off  [lioj  •  Errors  attendant  to 
determining  the  direction  cosines  of  this  direction  necessitate  degrading 
the  estimate  of  absolute  error  of  the  tantalum  results  to  +  1.0  percent. 
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MEASURED  ELASTIC  MODULI  [x  I012  dyne -cm’ 


2 

Fig.  1  Measured  values  of  ^>V  for  tantalum  as  a  function  of  temperature. 

The  direction  of  elastic  wave  propagation  for  this  sample  was  8° 
removed  from  [ll^  (See  Text). 
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TABLE  II.  The  elastic  constants  of  tantalum,  in  units  of  10  dyne-cm  , 
as  a  function  of  absolute  temperature.  Calculated  by  a  perturbation  method 
from  smooth  curves  drawn  through  values  obtained  from  measurements  made  in 
a  direction  of  propagation  8°  from  Q.10]  .  The  number  of  significant 
figures  shows  the  internal  consistency  of  the  data  and  does  not  indicate 
the  accuracy  of  the  absolute  value. 


<°K) 

cn 

C12 

C44 

2^Cll“C12^ 

3(C11+2C12> 

C12~C< 

0 

2.6632 

1.5816 

0.8736 

0.5408 

1.9421 

0.7080 

20 

2.6633 

1.5823 

0.8728 

0.5405 

1.9426 

0.7095 

40 

2.6641 

1 . 5844 

0.8709 

0.5398 

1.9443 

0.7135 

60 

2.6649 

1.5870 

0.8682 

0.5390 

1.9463 

0.7188 

80 

2.6652 

1.5896 

0.8647 

0.5378 

1.9481 

0.7249 

100 

2.6638 

1.5912 

0.8609 

0.5363 

1.9487 

0.7303 

120 

2.6603 

1.5909 

0.8564 

0.5347 

1  9474 

0.7345 

140 

2.6562 

1.5900 

0.8522 

0.5331 

1.9454 

0.7378 

160 

2.6518 

1.5888 

0.8480 

0.5315 

1.9431 

0.7408 

180 

2.6467 

1.5875 

0.8436 

0.5296 

1.9406 

0.7439 

200 

2.6412 

1.5860 

0.8393 

0.5276 

1.9377 

0.7467 

220 

2.6350 

1.5840 

0.8352 

0.5255 

1.9343 

0.7488 

240 

2.6286 

1.5818 

0.8309 

0.5234 

1.9307 

0.7509 

260 

2.6221 

1.5794 

0.8266 

0.5214 

1.9270 

0.7528 

280 

2.6155 

1.5768 

0.8224 

0.5194 

1.9230 

0.7544 

300 

2.6091 

1.5743 

0.8182 

0.5174 

1.9192 

0.7561 
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The  elastic  constant  o‘(C11+C«2)  +C, ,  for  tungsten  as  a  function  of  temperature 
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Fig.  3  The  sheer  constants,  C, ,  and  — (C. . -C.  „)  ,  for  tungsten  as  functions  of  temperature. 
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TABLE  III.  The  elastic  constants  of  tungsten,  in  units  of  10  dyne-cm  , 
as  a  function  of  absolute  temperature.  The  first  three  columns  are  readings 
taken  from  smooth  curves  drawn  through  the  experimentally  determined  points. 
The  number  of  significant  figures  shows  the  internal  consistency  of  the  data 
and  does  not  indicate  the  absolute  value. 


T 

(°K) 

evi2 

CL 

fV/ 

C44 

i  K 

2<circi2> 

cn 

C12 

B 

!<CU+2C12> 

C12'C44 

0 

5.3188 

1.6313 

1.6380 

5.3255 

2.0495 

3.1415 

0.4182 

20 

5.3188 

1.6313 

1.6380 

5.3255 

2.0495 

3.1415 

0.4182 

40 

5.3187 

1.6313 

1.6380 

5.3254 

2.0494 

3.1414 

0.4181 

60 

5.3187 

1.6312 

1.6379 

5.3254 

2.0496 

3.1415 

0.4184 

80 

5.3186 

1.6312 

1.6378 

5.3252 

2.0496 

3.1415 

0.4184 

100 

5.3181 

1.6310 

1.6368 

5.3239 

2.0503 

3.1415 

0.4193 

120 

5.3166 

1.6298 

1.6352 

5.3220 

2.0516 

3.1417 

0.4218 

140 

5  c  3138 

1.6273 

1.6315 

5.3180 

2.0550 

3.1427 

0.4277 

160 

5.3100 

1.6238 

1.6270 

5.3132 

2.0592 

3.1439 

0.4354 

180 

5.3036 

1.6207 

1.6220 

5.3049 

2.0609 

3.1422 

0.4402 

200 

5.2942 

1.6178 

1.6172 

5.2936 

2.0592 

3.1373 

0.4414 

220 

5.2846 

1.6153 

1.6124 

5.2817 

2.0569 

3.1318 

0.4413 

240 

5.2750 

1.6130 

1.6077 

5.2697 

2.0543 

3.1261 

0.4406 

260 

5.2653 

1.6108 

1.6031 

5.2576 

2.0514 

3.1201 

0.4406 

280 

5.2557 

1.6089 

1.5984 

5.2452 

2.0484 

3.1140 

0.4395 

300 

5.2462 

1.6072 

1.5937 

5.2327 

2.0453 

3.1078 

0.4381 
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Fig.  4  The  elastic  constant  — (C^j+Cj^)  +C44  for  "olybdenuni  as  a  function  °f  temperature.  Two  separate 
measurements  of  the  constant,  as  indicated  by  open  and  closed  circles,  were  made  between  liquid 
air  and  room  temperatures . 
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Fig.  5  The  shear  constants,  C,,  and  ~ci  *  for  ™olyb<ienu:m  *s  functions  of  temperature. 
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TABLE  IV.  The  elastic  constants  of  molybdenum,  in  units  of  10  dyne-cm  , 
as  a  function  of  absolute  temperature.  The  first  three  columns  are  readings 
taken  from  smooth  curves  drawn  through  the  experimentally  determined  points. 
The  number  of  significant  figures  shows  the  internal  consistency  of  the  data 
and  does  not  indicate  the  absolute  value. 


T 

(°K) 

fvi2 

CL 

£>V22 

C44 

ev32 

i(circi2> 

cu 

C12 

B 

i(Cu+2Ci2> 

C12~C44 

0 

4.3650 

1.2503 

1.3855 

4.5002 

1.7292 

2.6529 

0.4789 

20 

4.3639 

1.2497 

1.3855 

4.4997 

1.7287 

2.6524 

0.4790 

40 

4.3600 

1.2488 

1.3850 

4.4962 

1.7262 

2.6495 

0.4774 

60 

4.3557 

1.2479 

1.3840 

4.4918 

1.7238 

2.6465 

0.4759 

80 

4.3548 

1.2468 

1.3820 

4.4900 

1.7260 

2.6473 

0.4792 

100 

4.3495 

1.2456 

1.3807 

4.4846 

1.7232 

2 „ 6437 

0.4776 

120 

4.3389 

1.2446 

1.3770 

4.4713 

1.7173 

2.6353 

0.4727 

140 

4.3292 

1.2429 

1.3719 

4.4582 

1.7144 

2.6290 

0.4715 

160 

4.3212 

1.2402 

1.3672 

4.4482 

1.7138 

2 . 6253 

0.4736 

180 

4.3160 

1.2368 

1.3662 

4.4454 

1.7130 

2.6238 

0.4762 

200 

4  . 3146 

1.2333 

1.3638 

4.4451 

1.7175 

2.6267 

0.4842 

220 

4.3127 

1.2296 

1.3595 

4.4426 

1.7236 

2.6299 

0.4940 

240 

4.3061 

1.2257 

1.3540 

4.4344 

1.7264 

2.6291 

0.5007 

260 

4.2988 

1 . 2220 

1.3487 

4.4255 

1.7281 

2.6272 

0.5061 

280 

4.2905 

1.2188 

1.3444 

4.4161 

1.7273 

2.6236 

0.5085 

300 

4.2825 

1.2165 

1.3417 

4.4077 

1.7243 

2.6188 

0.5078 
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DISCUSSION  OF  RESULTS 


The  elastic  behavior  of  a  cubic  crystal  is  completely  specified  by 
the  two  shear  moduli,  C  =  and  C *  —  "|(Cxi“Ci2^  and  the  reciProcal  com“ 
pressibility  or  bulk  modulus,  B  =  .  c  is  a  measure  of  the  re¬ 

sistance  to  shear  deformation  across  the  (100)  plane  in  the  [010j  direct¬ 
ion  and  C1  is  the  resistance  to  shear  deformation  across  the  (110)  plane 
in  the  [lip]  direction.  B  represents  the  resistance  to  deformation  by  a 
uniform  hydrostatic  pressure. 

For  tantalum  the  computed  values  of  C  and  C'  are  shown  plotted  against 
temperature  as  the  solid  lines  of  Figs.  6  and  7.  Fig.  8  shows  a  corre¬ 
sponding  solid  line  for  B.  These  plots  for  C  and  C'  show  that  these 
quantities  vary  with  temperature  in  a  manner  similar  to  that  observed  for 
other  metals.  They  approach  0°K  with  zero  slope  and  depend  linearly  on 
temperature  for  a  good  portion  of  the  temperature  range.  However  B  as  a 

function  of  temperature  (Fig.  8)  shows  a  departure  from  previously  ob- 

o 

served  behavior  by  displaying  a  positive  slope  up  to  about  100  K,  then  a 
slope  reversal  and  negative  slope  up  to  room  temperature. 

In  obtaining  the  elastic  constants  from  a  perturbation  calculation  it 
would  be  expected  that,  if  the  anisotropy  is  not  a  strong  function  of  tem¬ 
perature,  the  perturbation  corrections  would  be  small  and  of  the  same  sign 
at  all  temperatures.  Thus  if  the  experimental  points  decrease  monotonically 
with  increasing  temperature  as  in  Fig.  1,  the  elastic  constants  which  only 
differ  by  small  amounts  from  the  experimental  data  should  exhibit  the  same 
temperature  dependence.  That  this  is  not  the  case  for  B  is  shown  by  the 
solid  line  in  Fig.  8  which  is  plotted  from  the  elastic  constants  listed  in 
Table  II.  It  is  possible  that  a  fortuitous  choice  of  propagation  direction 
combined  with  the  inevitable  errors  in  angle  and  transit  times  may  lead  to 
the  observed  behavior. 

In  order  to  observe  the  effects  of  errors  of  measurement  or  orientation 
on  the  approximation  method  used  in  obtaining  the  elastic  constants  of 
tantalum,  the  input  data  were  slightly  changed  and  the  elastic  constants 

2 

recalculated  using  a  CDC  1604  computer.  Thus,  the  measured  values  of  ^>V 
were  individually  and  collectively  changed  by  +  1%,  the  polar  angles  being 
held  constant,  or  a  measured  polar  angle  0  or  0  was  changed  by  +  2°,  the 
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ig.  6  The  calculated  shear  constant  C, ^  for  tantalum  as  a  function  of  tem¬ 
perature  (solid  line).  The  dashed  lines  show  versus  T  when  (>V2, 
0,  or  0  are  changed  by  the  indicated  amounts. 


Fig.  7  The  calculated  shear  constant  ^or  tanta^uin  as  a  function  of 

temperature  (solid  line)  .  The  dashed  lines  show  ^(C^-C^)  versus  T 
when  ^>V2,  0,  or  0  are  changed  hy  the  indicated  amounts. 


ig.  8  The  calculated  bulk  modulus,  j(C^+2C^2)»  ^or  tantalura  as  a  function 

of  temperature  (solid  line).  The  dashed  lines  show  versus 

T  when  (pV^,  0,  or  0  are  changed  by  the  indicated  amounts. 


values  of  ^>V  and  the  other  angle  being  held  constant.  The  results  of 
these  calculations  are  indicated  as  dashed  lines  in  Figs.  6,  7,  and  8 
and  show  that  for  substantial  input  data  changes  the  curves  for  C,  C1 , 
and  B  are  merely  shifted  up  or  down  their  respective  ordinates  with  no 
changes  in  shape.  Because  there  is  no  change  of  shape  with  these  per¬ 
turbations  of  the  input  data,  it  is  concluded  that  the  observed  slope 
reversal  effect  in  the  bulk  modulus  curve  for  tantalum  is  a  real  effect. 

In  the  tungsten  crystal  C  and  C'  were  measured  directly  since  pro¬ 
pagation  was  along  [lio]  .  These  results  are  shown  in  Fig.  3.  The 
variation  of  B  with  temperature  has  been  calculated  and  is  shown  as  the 
upper  curve  in  Fig.  9. 

The  most  striking  feature  of  the  variations  of  the  measured  quanti¬ 
ties  C,  C'  and  2^CH+C12^  +C44  temPerature  is  their  very  flat  be¬ 

havior  from  liquid  air  down  to  liquid  helium  temperatures.  This  is  the 
result  of  there  being  no  observable  difference  in  the  elastic  wave  transit 
times  for  the  three  modes  from  liquid  air  to  liquid  helium  temperatures. 
Because  of  this,  the  bulk  modulus,  B,  remains  constant  with  increasing 
temperature  out  to  approximately  100°K.  Beyond  this  point  B  displays  a 
positive  slope  out  to  160°K  when  it  takes  up  a  negative  slope  and  a  linear 
variation  with  temperature. 

Fig.  3  shows  tungsten  is  nearly  isotropic  in  the  temperature  region 
measured,  being  exactly  so  at  the  cross-over  point  between  the  two  modes 
at  about  195°K.  Above  140°K  the  shear  modes  vary  quite  linearly  with 
temperature . 

As  was  the  case  with  tungsten,  C  and  C1  were  measured  directly  for 
molybdenum.  These  results  are  displayed  in  Fig.  5,  which  shows  that  C 
varies  monotonically  with  temperature  but  C1  displays  a  tendency  to  level 
off  in  the  90°K  and  180°K  regions. 

The  variation  of  B  with  temperature  has  been  calculated  and  is  shown 
as  the  lower  curve  in  Fig.  9.  It  is  seen  that  the  slope  of  the  bulk  modulus 
curve  goes  from  zero  to  negative  as  the  temperature  rises  from  0°K  but  there 
are  two  subsequent  pairs  of  gradual  slope  reversals  over  the  remainder  of 
the  temperature  range  measured,  one  reversal  at  about  80°K  and  one  at  about 
230°K.  Reference  to  the  1%  error  flags  in  Figs.  5  and  9  readily  establishes 
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Fig.  9  The  adiabatic  bulk  moduli  tCCji+ZCjo) >  for  tungsten  and  molybdenum  as  functions  of  temperature. 


that  these  slope  variations  of  C'  and  B  all  lie  within  about  0.2%  of  a 
line  that  would  show  monotonically  decreasing  behavior  with  temperature. 
However  as  Fig.  4  indicates,  two  separate  measurements  of  the  combination 
+C44  versus  temperature  (the  longitudinal  acoustic  mode  for 
L110J  )  were  made  from  liquid  air  to  room  temperatures,  with  consistent 
results.  This  consistency  indicates  that  the  depicted  rather  than  the 
monotonic  behavior  of  B  versus  T  is  warranted. 

It  is  seen  from  Fig.  1-5  that  the  temperature  variation  of  the  ob¬ 
served  quantities  is  similar  to  that  observed  for  other  metals.  That  is 
the  elastic  constants  approach  0°K  with  zero  slope  and,  except  for  molybdenum, 
are  remarkably  linear  over  a  wide  temperature  range.  The  percentage  change 
of  the  constants  relative  to  their  0°K  values  is,  with  the  exception  of 
tantalum,  about  the  same  as  that  observed  for  the  noble  metals.  Near  room 
temperature  all  the  constants  vary  linearly  with  temperature  and  have 
significantly  smaller  logarithmic  derivatives  than  the  noble  metals.  Un¬ 
fortunately,  the  change  of  elastic  constants  with  pressure  have  not  yet 
been  measured  for  these  materials  so  that  it  is  not  possible  to  evaluate 

the  intrinsic  temperature  dependence  of  the  elastic  constants. 

8  9 

Other  investigators  *  have  reported  results  for  these  metals  obtained 
solely  by  dynamic  methods.  Also  reported  in  the  literature^*  ^  are  re¬ 
sults  obtained  by  dynamic  measurements  methods  combined  with  Bridgman's 
12 

compressibility  data. 

Table  V  compares  room  temperature  values  from  the  literature  with  the 
results  of  this  paper.  The  differences  between  the  tantalum  results  re¬ 
ported  by  Bolef  and  by  this  paper  tend  to  be  close  to  the  sum  of  the  re¬ 
spective  estimates  of  error  ( '—  2%) .  This  represents  fair  agreement.  The 
tungsten  results  reported  by  Bolef  and  de  Klerk  are  in  quite  good  agree¬ 
ment  with  those  reported  in  this  paper  (within  1%) . 

We  note  that  for  tungsten  the  difference  between  Bridgman's  value  of 
the  compressibility  and  the  result  of  this  work  is  very  close  to  be  cal¬ 
culated  difference  between  the  room  temperature  adiabatic  and  isothermal 
compressibilities.  For  tantalum  and  molybdenum,  however,  no  such  agree¬ 
ment  is  found,  using  either  set  of  ultrasonically  derived  constants. 
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TABLE  V.  The  reported  room  temperature  values  of  the  adiabatic  elastic 
constants  of  tantalum,  tungsten,  and  molybdenum,  in  units  of  10^  dyne-cm 


Element 

C 

C' 

B 

Investigator 

Tantalum 

0.825 

0  o  53 

1.96 

Bolef3 

0.818 

0.518 

1  919 

This  Paper 

Tungsten 

1.51 

1.51 

2.99 

Wright^ 

1.53 

1.54 

3.08 

Bridgman0 

1.604 

1.598 

3.083 

Bolef  &  de  Klerkd 

1.607 

1,594 

3.108 

This  Paper 

Molybdenum 

1.09 

1.4 

2.73 

e 

Druyvesteyn 

1.068 

1.510 

2.683 

Bolef  6c  de  Klerkd 

1.216 

1.342 

2.619 

This  Paper 

a 

See 

Reference 

8 

b 

See 

Reference 

11 

c 

See 

Reference 

12 

d 

See 

Reference 

9 

e 

See 

Reference 

10 
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The  difference  in  the  values  of  the  elastic  shear  constants  reported 

for  molybdenum  by  Bolef  and  de  Klerk  and  by  this  paper  are  quite  large. 

It  seems  clear,  particularly  in  light  of  the  agreement  on  tungsten,  that 
these  differences  in  moduli  are  real  and  that  they  depend  on  physical 
differences  between  the  specimens  rather  than  on  differences  in  the  mea- 

suring  apparatus.  That  this  must  be  the  case  is  made  more  credible  by 

considering  that  in  measuring  C^,  which  is  measured  directly  for  [no] 
propagation,  the  difference  in  the  two  reported  values  corresponds  to  a 
difference  in  transit  time  of  0.5  ^us  in  the  sample  used,  far  outside  the 
region  of  any  reading  error. 

A  possible  explanation  of  the  difference  in  reported  room  temperature 
elastic  constants  (~,12%  differences  for  C  and  C'  from  the  values  of  Bolef 
and  de  Klerk)  for  molybdenum  is  the  influence  of  alloying  of  the  impurities. 
Small  amounts  of  a  second  phase  have  been  found  to  change  the  single  crystal 
elastic  constants  of  copper , silver, ^  and  lithium"^.  These  changes 
in  elastic  constant  have  been  interpreted  in  terms  of  changes  in  the  long 
range  and  short  range  contributions  to  the  elastic  constants.  The  experi- 
ments  indicate  that  the  percent  change  in  modulus  per  percent  change  in 

composition  is  seldom  more  than  four  and  usually  near  unity.  Since  for 

molybdenum  the  observed  percent  change  in  modulus  per  percent  change  in 

composition  is  about  100,  it  is  difficult  to  see  how  these  differences  can 

be  interpreted  on  the  basis  of  simple  alloying. 

Anelastic  effects  due  to  dislocations  also  contribute  to  change  the 

moduli.  A  dislocation  contribution  (diminution)  to  the  moduli  of  about  one 

17  18 

percent  has  been  observed  for  a  single  crystal  of  copper  and  of  lead 
of  purity  in  the  10  to  200  ppm  range.  In  these  experiments,  the  addition 

of  small  amounts  of  impurity  tended  to  raise  the  observed  modulus  through 

19 

dislocation  pinning.  Chambers  and  Schultz  have  observed  an  internal 
friction  spectrum  for  molybdenum  in  the  low  temperature  range  that  shows  a 
pronounced  peak  in  the  region  just  below  100°K  (the  peak  shifts  with 
annealing)  and  a  high  level  minimum  at  about  200°K.  They  explain  certain 
of  the  observed  internal  friction  effects  as  being  due  to  the  interactions 
of  gaseous  interstitials  and  dislocations.  The  assay  of  the  present  moly¬ 
bdenum  sample  (Table  I)  shows  a  significantly  higher  gas  content  than  the 
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sample  of  Bolef  and  de  Klerk.  In  addition,  the  present  sample  is  much 

larger,  and  since  dislocation  densities  increase  rapidly  with  the  dia= 

20 

meter  of  the  specimen  ,  conditions  are  right  in  the  present  sample  for 
an  anelastic  mechanism  such  as  that  proposed  by  Chambers  and  Schultz  to 
be  quite  pronounced. 

The  differences  in  the  reported  elastic  shear  constants  for  moly¬ 
bdenum  are  in  opposite  directions,  and  if  these  differences  can  be  ex¬ 
plained  by  means  of  an  anelastic  mechanism,  then  these  differences  should 

yield  some  evidence  as  to  how  dislocations  are  distributed  in  the  lattice. 

21 

Seeger  has  conducted  an  analysis  of  low  temperature  internal  friction 

peaks  wherein  only  a  small  fraction  of  the  dislocations  are  considered  to 

22 

lie  along  crystallographic  directions.  The  work  of  Thompson  and  Holmes 
on  copper  indicates,  however,  that  an  appreciable  fraction  of  the  dislo¬ 
cations  tend  to  lie  along  principal  crystallographic  directions.  That 
the  two  orthogonal  shear  modes  in  the  molybdenum  sample  are  almost  dia¬ 
metrically  affected  by  this  presumed  dislocation-interstitial  gas  anelastic 
mechanism  indicates  that  such  a  preferred  alignment  of  dislocations  may 
indeed  be  present. 
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BULK  MODULUS 


As  shown  in  Fig.  8  the  bulk  modulus  of  tantalum  rises  with  increasing 
temperature  above  the  4.2°K  value  before  beginning  to  decrease  linearly. 
Observation  of  this  anomalous  behavior  of  B  with  temperature  has  led  the 
authors  to  calculate  the  bulk  moduli  of  the  other  transition  elements 

whose  elastic  constants  have  been  determined  ultrasonically  over  the  same 

23 

temperature  range.  The  following  data  were  used:  Alers  for  vanadium; 

24  25 

Rayne  and  Chandrasekhar  for  iron;  Alers,  Neighbours,  and  Sato  for 

26  27 

nickel;  Rayne  for  palladium;  Carroll  for  niobium;  and  this  paper  for 
tantalum,  tungsten,  and  molybdenum.  Curves  presenting  these  results  are 
shown  in  Figs.  10  and  11  where  the  bulk  modulus  has  been  normalized  to  its 
0°K  value,  denoted  by  Bq.  Fig.  10  shows  B/Bq  versus  temperature  for  the 
body-centered  cubic  transition  metals.  Here  the  behavior  of  niobium, 
tantalum,  tungsten,  and  molybdenum  is  significantly  different  from  that 
displayed  by  vanadium  and  iron.  Similarly  in  the  f ace-centered  cubic 

elements  shown  in  Fig.  11  palladium  has  quite  a  different  temperature  de~ 

26 

pendence  from  that  displayed  by  nickel.  Rayne  has  ascribed  the  departure 
from  linearity  at  temperatures  above  100°K  of  the  elastic  constants  of 
palladium  to  the  change  with  temperature  of  the  contribution  to  the  elastic 

constants  from  the  unfilled  d  band.  B/B  values  from  the  elastic  con- 
28  1  ° 

stants  of  copper  ,  silver  and  gold  have  been  calculated  and  are  included 
in  Fig.  11  to  indicate  the  temperature  dependence  of  the  bulk  moduli  of 
metals  that  have  closed-d  shells. 

Taking  monotonically  decreasing  behavior  with  temperature  as  the  norm, 
Figs.  10  and  11  show  that  the  3d  and  closed-d  shell  metals  follow  the  norm. 
Niobium  and  palladium,  4d  metals,  change  the  sign  of  their  slope  several 
times  over  the  temperature  range  measured,  displaying  a  strong  positive 
slope  as  room  temperature  is  approached.  Tantalum  (5d)  shows  an  early 
positive  slope,  then  tends  to  follow  the  norm  above  100°K.  Tungsten  (5d) 
has  a  flat  response  from  0°K  out  to  100°K  where  a  positive  slope  finally 
asserts  itself.  B/Bq  for  tungsten  then  parallels  that  of  tantalum  very 
closely  on  up  to  room  temperature.  Molybdenum  (4d)  generally  follows  the 
norm  but  with  undulating  departures. 

Figs.  10  and  11  indicate  a  significant  difference  between  the  behavior 
of  the  bulk  moduli  of  the  3d  metals,  and  the  4d-5d  metals.  Such  a  dif- 
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Pig.  10  The  bulk  nodull,  normalized  to  their  respective  0  K  values,  for  certain  of  the  body-centered 
cubic  transition  metals  as  functions  of  temperature.  Solid  lines  are  used  for  3d  metals, 
dashed  lines  for  4d  and  3d  metals. 
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Fig.  11  The  bulk  moduli,  normalized  to  their  respective  0°K  values,  for  certain 
of  the  face-centered  cubic  transition  metals  and  for  the  noble  metals. 


ference  is  in  consonance  with  the  observance  of  ferromagnetism  in  the  3d 
group  but  not  the  4d=5d  groups. 

One  possible  reason  for  this  anomalous  behavior  of  B  with  temperature 
that  suggests  itself  is  that  these  metals  are  involved  in  antiferromagnetic 
ordering  in  these  temperature  ranges,  with  associated  anomalies  in  their 

thermal  expansion  coefficients.  Evidence  in  support  of_  this  from  other  ex- 
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periments  is  scarce  and  difficult  to  assess.  Bolef,  de  Klerk,  and  Brandt 

have  reported  anomalous  thermal  expansion  behavior  in  the  180°K  to  240°K 

region  for  vanadium.  Fig.  10,  however,  shows  that  vanadium  is  a  well-behaved 

30 

3d  metal  in  its  B  behavior.  Hoare  and  Matthews  found  a  maximum  in  the 
magnetic  susceptibility  of  palladium  between  80°K  and  100°K  yet  Crangle  and 
Smith  recently  looked  for  a  specific  heat  A  -point  in  this  region  for 

palladium,  with  no  success.  They  also  measured  the  specific  heat  of  vana- 
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dium,  confirming  the  results  of  Clusius,  Franzosini,  and  Piesbergen  that 
vanadium  shows  no  anomaly  between  175°K  and  265°K,  the  region  of  Bolef, 
el  al's  anomalous  thermal  expansion  behavior  for  the  same  metal.  From  the 
point  of  view  that  the  body-centered  cubic  structure  is  the  more  natural 
vehicle  for  exhibiting  antiferromagnetic  behavior,  and  looking  at  the 
anomalous  behavior  of  B/Bq  versus  T  for  the  4d-5d  members  of  this  group 
as  possible  evidence  of  this  effect,  then  it  is  somewhat  disconcerting  that 
face-centered  cubic  palladium  displays  the  same  anomalous  behavior.  As  the 
results  of  detailed  specific  heat,  magnetic  susceptibility,  and  thermal  ex¬ 
pansion  experiments  for  pure  samples  of  these  metals  become  available,  the 
role  of  antiferromagnetic  ordering  in  the  transition  metals  may  become 
clearer . 
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For  body-centered  and  f ace°centered  cubic  structures  de  Launay  has 
raised  the  possibility  that  is  equal  to  the  bulk  modulus  of  the 

electron  gas  of  a  metal.  The  temperature  dependence  of  a  free  electron 
bulk  modulus  should  show  a  steady  decrease  with  increasing  temperature  as 
the  lattice  parameter  increases.  Reference  to  the  C12~C44  co^umns  in 
Tables  II,  III,  and  IV  shows  that  tantalum,  tungsten,  and  molybdenum  like 
copper  and  silver  exhibit  just  the  opposite  effect,  an  increase  of  c-j.2”C44 
with  temperature.  This  reinforces  de  Launay' s  original  admonition  that  the 
bulk  modulus  as  derived  from  his  considerations  is  a  phenomenological  quantity. 
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It  does  not  reflect  on  the  validity  of  the  Sommerfeld  free  electron  model 

However  correlation  of  C,„-C. ,  with  the  electronic  bulk  modulus  of  mono- 

12  44  34 

valent  metals  is  still  considered  a  possibility 

Excluded  in  the  above  calculations  is  the  energy  of  the  lowest  state 

Wq.  Although  its  effect  on  the  bulk  modulus  is  usually  taken  to  be  small 

it  is  possible  that  this  may  not  be  the  case.  Further  consideration  of 

this  contribution  await  a  cellular  method  calculation  of  W  , 

o 
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OVERLAP  FORCES 


35 

Fuchs  has  accounted  for  the  elastic  constants  of  copper  as  being 
composed  of  two  main  parts:  the  electrostatic  stiffnesses  due  to  that 
part  of  the  energy  which  can  be  regarded  as  the  energy  of  an  array  of 
positive  point  charges  immersed  in  a  uniform  negative  sea;  and  the  over¬ 
lap  energy  due  to  the  overlap  of  the  charge  distributions  of  ions  plus 

36  37 

the  exchange  energy.  Following  the  lead  of  Fuchs  and  Zener  ,  Isenberg 
has  concluded  that  the  major  contribution  to  the  elastic  constants  of  the 
body-centered  cubic  transition  metals  is  due  to  the  overlap  energy  (the 
Coulomb  interaction  energy  due  to  the  overlap  of  the  charge  distribution 
of  ions  plus  the  exchange  energy  due  to  the  same  overlap). 

The  relations  giving  the  overlap  contributions  to  the  elastic  con- 

g 

stants  can  be  written  as: 

N  t  V]  i  =  3/8  (C44-C12) 

N  tV  3  1  =  3/4  (C12+2C44> 

N  tV]  2  =  CirC44 

where  N  is  the  density  of  atoms,  is  the  nearest  neighbor  interaction 
energy  and  is  the  next  nearest  neighbor  interaction  energy,  and  the 

primes  represent  differentiation.  To  a  first  approximation,  the  first 

37 

derivative  of  may  be  neglected  ,  and  the  remaining  derivatives  are 
evaluated  at  the  nearest  of*  next=nearest  neighbor  distance  as  appropriate. 

In  applying  the  above  formulas,  only  the  part  of  the  elastic  constants  due 
to  overlap  forces  is  included,  i.e.,  the  overlap  contribution  to  the  modu¬ 
lus  is  taken  to  be  the  observed  modulus  minus  the  contributions  from  all 
other  sources. 

Using  the  room  temperature  values  of  the  elastic  constants  of  tungsten 
as  measured  by  Wright11,  Isenberg  determined  that  the  second  derivative  of 
the  next-nearest  neighbor  interaction,  W^' 1 ,  is  positive,  indicating  a  re¬ 
pulsive  second-nearest  neighbor  interaction  in  tungsten.  The  values  he 
obtained  for  the  first  and  second  derivatives  of  the  nearest  neighbor  in¬ 
teraction  W1’  and  W^"  indicated  a  binding  interaction  for  nearest  neighbors. 


29 


Isenberg  interpreted  these  results  as  a  possible  substantiation  of  Zener's 
3  6 

proposal  of  an  antiferromagnetic  arrangement  for  the  body-centered  cubic 
transition  elements  vanadium,  chromium,  niobium,  molybdenum,  tantalum,  and 
tungsten. 


Taking  the  published  values  of  the  0  K  elastic  constants  of  vanadium  , 
27  24 

niobium  ,  and  iron  and  the  results  of  this  paper  these  same  derivatives 

11  -2 

have  been  calculated.  The  results  in  units  of  10  dynes  cm  are  shown 


in  Table  VI .  For  all  these  transition  elements  except  iron  the  results  are 
similar  to  Isenberg' s.  That  is,  the  first  derivative  of  the  nearest  neighbor 
interaction  energy  is  negative  and  the  second  derivatives  of  both  interaction 
energies  are  positive. 

The  magnitudes  of  the  second  derivatives  generally  increase  as  one 
passes  down  and  to  the  right  in  the  periodic  table.  The  magnitude  of  the 
first  derivative  of  increases  to  the  right  in  the  periodic  table.  The 
derivatives  for  iron  are  different  in  sign  and  magnitude  as  might  be  ex¬ 
pected  for  a  ferromagnetic  material. 

°OL  r 

We  note  that  if  a  simple  exponential  form  -  A^e  i  is  chosen  for 
both  and  W2  then  without  neglecting  W^'  it  can  be  shown 


where  the  subscript  on  the  square  bracket  denotes  the  appropriate  distance 
at  which  the  quantity  in  the  bracket  is  to  be  evaluated.  The  elastic  con¬ 
stants  used  above  are  understood  to  include  only  that  part  due  to  overlap. 
For  the  transition  metals  considered  here  anc*  ^12^^  so  t^iat  t^ie 

above  relations  predict  W^  and  W2  both  to  be  repulsive  energies.  This  re¬ 
sult  is  in  agreement  with  ascribing  the  large  cohesive  energy  of  molybdenum 

38  39 

and  tungsten  to  the  large  number  of  electrons  in  the  low  lying  d  band. 
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TABLE  VI.  Values  of  the  derivatives  of  the  nearest  and  next-nearest 
neighbor  interaction  energies  attributable  to  overlap  forces,  as 
calculated  from  the  0°K  elastic  constants.  Here  N  =  the  density  of 
atoms,  =  the  nearest  neighbor  interaction  energy,  =  the  next- 
nearest  neighbor  interaction  energy  and  the  primes  indicate  differ¬ 
entiation  with  respect  to  r.  The  subscripts  1  and  2  outside  the 
brackets  indicate  evaluation  of  the  enclosed  quantities  at  the 

nearest  and  next-nearest  neighbor  distances,  respectively.  The  value 
11  -2 

given  are  in  units  of  10  dyne-cm 


Element 

N  [r  v] 

1 

N  |jr2W1 '  ' 

1. 

N  [r2W2" 

1. 

Vanadium 

-1.8 

9.2 

15.4 

Niobium 

-3.2 

10.2 

20.2 

Tantalum 

-2.1 

20.5 

15.9 

Molybdenum 

-1.0 

26.2 

29.9 

Tungsten 

-0.8 

34.5 

34.5 

Iron 

+0.8 

20.2 

7.8 
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DEBYE  TEMPERATURES 


Using  Houston's  method  as  outlined  by  Betts,  Bhatia,  and  Wyman^  the 
Debye  temperatures  of  tantalum,  tungsten,  and  molybdenum  have  been  calcu¬ 
lated  from  their  0°K  elastic  constants.  Table  VII  presents  these  results 
and  compares  them  with  the  results  of  calorimetric  measurements.  Although 
the  values  entered  for  the  Debye  temperatures  of  tantalum  in  Table  VII 
differ  by  rather  large  amounts,  the  lattice  specific  heat  given  by 

G  =  261. 9°K  is  not  in  disagreement  with  the  lower  temperature  part  of  the 

41  42 

specific  heat  data  of  White  et  al  and  Worley  et  al. 

43 

The  data  of  Keesom  and  Desirant  are  too  scattered  to  admit  critical 
44  o 

comparison  and  Wolcott  lists  only  the  value  9q  =  245  K.  Thus  the 

situation  for  comparison  of  elastic  and  calorimetric  Debye  temperatures  of 

45 

tantalum  is  similar  to  that  for  other  metals  ;  i.e„,  where  the  calorimetric 
data  is  not  too  scattered,  the  lattice  specific  heat  calculated  from  the 
elastic  constants  fits  the  low  temperature  part  of  the  calorimetric  data. 

Much  the  same  remarks  apply  to  tungsten  and  molybdenum.  For  tungsten 

46 

the  data  of  Clusius  and  Franzosini  is  fitted  at  the  low  temperature  end 

by  G  s  384. 3°K.  Clusius  and  Fronzosini  list  G  =  380 °K  and  both  Debye 

°  47  ° 

temperatures  fit  the  data  quite  well.  As  in  the  case  for  tantalum  Wolcott 

lists  only  a  number  and  other  data  in  the  literature  is  too  scattered  for 

comparison. 

For  molybdenum,  the  calculated  value  Gq  =  47402°K  gives  a  lattice 
specific  heat  in  agreement  with  the  lower  temperature  portion  of  the  data 
of  Clusius  and  Franzonini.  The  other  published  calorimetric  data  is  not 
suitable  for  comparison. 
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TABLE  VII.  The  Debye  temperatures  of  tantalum,  tungsten,  and  molybdenum 
as  calculated  from  the  elastic  constants  and  as  derived  from  low  temperature 
specific  heat  data. 


Element 

9 

o 

Reference 

Tantalum 

261.9 

This  Work 

255 

a 

246.5 

b 

245 

c 

231.0 

d 

Tungsten 

384  o  3 

This  Work 

380 

e 

405 

c 

Molybdenum 

474  o  2 

This  Work 

470 

c 

454 

e 

445 

f 

a  See  Reference  41 
b  See  Reference  43 
c  See  Reference  44 
d  See  Reference  42 
e  See  Reference  46 

f  M.  Horowitz  and  J.G.  Daunt,  Phys  Rev  91,  1099  (1953) 
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